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Abstract 

Accurate representation of the physical layer is required for analysis and simulation of multi-hop networking in sensor, ad hoc, 
and mesh networks. This paper investigates, models, and analyzes the correlations that exist in shadow fading between links in 
multi-hop networks. Radio links that are geographically proximate often experience similar environmental shadowing effects and 
thus have correlated fading. We describe a measurement procedure and campaign to measure a large number of multi-hop networks 
in an ensemble of environments. The measurements show statistically significant correlations among shadowing experienced on 
different links in the network, with correlation coefficients up to 0.33. We propose a statistical model for the shadowing correlation 
between link pairs which shows strong agreement with the measurements, and we compare the new model with an existing shadowing 
correlation model of Gudmundson (1991). Finally, we analyze multi-hop paths in three and four node networks using both correlated 
and independent shadowing models and show that independent shadowing models can underestimate the probability of route failure 
by a factor of two or greater. 

Index Terms 

Wireless sensor, ad hoc, mesh networks, shadowing, correlation, statistical channel model, wireless communication, measure- 
ment, performance 

I. Introduction 

Both simulation and analysis are critical to the development of multi-hop networks, including mesh, ad hoc, and sensor net- 
works. However, current physical layer models do not accurately represent radio channels in multi-hop wireless networks HI, 
and as a result, there is a significant disconnect between simulation and real world deployment. There is significant interest in im- 
proving statistical models beyond the current state-of-the-art in order to decrease the difference between simulation and analysis 
results and experimental deployment results. 

This paper presents a statistical joint path loss model between a set of static nodes. Joint path losses and transmit powers 
determine the connectivity, reliability in interference, and energy consumption during power control, of network communications. 
Channel models used in multi-hop networks have considered path losses to be independent, yet they are correlated through 
shadowing effects. We demonstrate these correlations via measurements and present a correlated shadowing loss model, which is 
then shown to have a dramatic effect on network connectivity. 

We do not address other random processes like transmit power variation, manufacturing variations among nodes, position 
of nodes in random deployments, mobility of nodes, or interference models. However, the developed model informs future 
development of path loss models for mobile networks, and may be used to analyze the effects of other variation and interference 
models. 

A. Single-Link Path Loss Model 

Radio propagation measurement and modeling for a single radio link has been reported extensively over the past century 121, 
13, m. 111. In general, when there is no site-specific knowledge of the environment, the ensemble mean received power, P{d) 
(dBm), at a distance d from the transmitter, is [Si, H, 

P{d) = Pt - Ho - lOnp logio X-' (1) 

where Pt is the transmitted power in dBm, Up is the path loss exponent, and Ho is the loss experienced at a short reference 
distance Aq from the transmitter antenna. This model incorporates the free space path loss model when rip = 2, and extends to 
practical (obstructed) multipath environments when rip > 2. 

On a particular link, received power will vary from the ensemble mean because of fading. The measured received power for the 
link between transmitter i and receiver j is, 

P,j ==P(d,,,)-Z,,„ (2) 



where di,j is the distance between nodes i and j, and Zij is the fading loss. In general, shadow fading, small-scale or frequency- 
selective fading, and antenna and device losses all contribute to Zij. Wideband receivers reduce the effects of small-scale or 
frequency-selective fading issues, and antenna and device-caused variations are generally small compared to shadowing variations. 
Shadow fading, also called medium-scale fading Q, describes the loss suffered as the signal passes through or diffracts around 
major obstructions in its path from the transmitter to the receiver These obstructions include walls and furniture indoors, and 
buildings, terrain, and trees outdoors. 

We hypothesize that shadowing losses are correlated across different links which are geographically proximate. Since shadow- 
ing is central to the analysis in this paper, we separate total fading loss Zij into two contributions, 

Zij = Xij + Yi_j , (3) 

where Xij represents the shadowing loss, and Yij represents all other (non-shadowing) losses. 

B. Application in Multi-hop Networking Research 

In the multi-hop networking simulation and analysis literature, two models are used: 

1) The circular coverage model: Zij =- for all links, and thus the coverage area is a perfect circle, as shown in Figurellta). 

2) The i.i.d. log-normal shadowing model: For all links (z, j), random variables Zj.j are independent and identically distributed 
Gaussian with zero mean and variance a\, as shown in Figure[TJc). 

We argue that both models are at opposite extremes, and both problematic. Note that 'reaUstic coverage' is commonly depicted 
pictorially as a coverage area with random range as a function of angle ||6|, iQ, as in Figure [TJb), and neither fading model 
produces such a random shape. It is easy to recognize that the deterministic, circular coverage areas are unrealistic for wireless 
communications links. However, circular coverage has been a common assumption in ad hoc and sensor network research and 
has been used to generate foundational research results. Kotz, Newport, and Elliot fS] examined the set of papers in the MobiCom 
proceedings from 1995 through 2002, and found that out of 36 papers which required radio models, only four did not use a circular 
coverage model. 

In comparison, the i.i.d. shadowing model is non-deterministic, and eliminates the concept of coverage area. Since the model 
has no spatial memory, even two nearly overlapping links would be represented as statistically independent. For example node 2 
in Figure [TJc) may be connected while node 1 is not. 

Recent research, including Hekmat and Van Mieghem |f7| and Bettstetter and Hartmann ||6l, has studied connectivity in ad 
hoc networks using the i.i.d. log-normal shadowing model. Their analyses indicate that for a constant level of connectivity, node 
deployment density can be reduced when the variance of the shadowing is increased. This increase in connectivity is largely a 
result of the model's independence assumption. Since losses in links in the same direction from a transmitter are independent, if 
one link is disconnected because of high loss, another node in the same direction is likely to be connected. 

In reality, if an obstacle in one direction from a transmitter strongly attenuates its signal, any receiver behind the obstacle is 
likely to experience high fading loss. For example, if the environment in Figure |2] causes severe shadowing, it is likely to cause 
additional path loss on both links a and h. In contrast, the i.i.d. log-normal shadowing model assumes that the shadowing across 
links a and 6 will be independent and thus exaggerate the connectivity. We quantify this argument in Section lVlIl 



C. Correlation Limits Link Diversity 

Diversity methods are common means to achieve reliability in unreliable channels. Multi-hop networking serves as a network- 
layer diversity scheme by allowing two nodes to be connected by any one of several multi-hop paths. All diversity schemes are 
limited by channel correlations. Correlations have been studied and shown to limit diversity gains in time, space, frequency and 
multipath diversity schemes ID, 13], H, ifTOll . 

Yet little research has addressed channel correlations on links in sensor, mesh, and ad hoc networks. This paper presents an 
initial investigation into quantitatively assessing the correlation in the shadow fading experienced on the different links of a multi- 
hop network. This investigation is experimental, using full link measurements of an ensemble of deployed networks to estimate 
and test for statistical correlations. We propose a joint path loss model which accurately represents observed correlations in link 
shadowing. Further, we quantify the effect that such correlation has on source to destination path statistics. We show that for a 
simple three node network that the probability of path failure is double what would be predicted by the i.i.d. log-normal shadowing 
model. 

II. Related Work 

Shadow fading correlations have been measured and shown to be significant in other wireless networks. For example: (1.) in 
digital broadcasting, links between multiple broadcast antennas to a single receiver have correlated shadowing which affects the 
coverage area and interference characteristics fTTJ : (2.) in indoor WLANs correlated shadowing is significant (as high as 0.95) 
strongly impacts system performance lfT2l : and (3.) in cellular radio correlation on links between a mobile station and multiple 
base stations significantly affects mobile hand-off probabilities and co-channel interference ratios lfT3l . lfT4l . ifTsl . 



In cellular radio, the model of Gudmundson fT6l is used to predict shadowing correlation for the link between a mobile station 
(MS) to a base station over time as the MS moves. In Section [VTl we address the difficulty in applying this model to multi-hop 
networks. We quantitatively compare it with the proposed model when the Gudmundson model may be applied. Wang, Tameh, 
and Nix ITTI extended Gudmundson's model to the case of simultaneous mobility of both ends of the link, for use in MANETs, 
and relate a sun-of-sinusoids method to generate realizations of the shadowing process in simulation. Both works use "correlated 
shadowing" to refer to the correlation of path loss in a single link over time, while the present work studies the correlation of many 
disparate links at a single time. 

The closest study to the present work used RSS measurements in a single network to quantify correlations between two links 
with a common node IfTSl . Those results could not be complete because a single measured network cannot provide information 
about an ensemble of network deployments. The present study uses multiple measured networks to examine many pairs of links 
with the identical geometry, both with and without a common node. 

III. Measurement Setup 

In this section we present our method for measuring the path-loss of each pair-wise link in a multi-hop network, using a 
specialized sensor network. This system is referred to as the network channel measurement system (NCMS). The NMCS allows 
us to quickly measure the received power Pij (dBm) of every link (z, j) in a deployed network, to measure across a range of 
frequencies, and to record the data on a laptop for later analysis. 

A. Equipment 

The nodes used in the measurement campaign are "mica2" motes manufactured by Crossbow. A mica2 mote operates in the 
902 - 928 MHz band using a Chipcon CCIOOO FSK transceiver. The transmit power is user programable and can be varied based 
on the network topology and environmental density. The mica2 measures and reports RSS values for each received signal |fT9l . 

1) Battery Variations: Transmit power is proportional to the battery voltage squared. Since measurement and data collection 
of one deployed network take on the order of minutes, battery voltage can be considered constant during each experiment. Each 
device measures and reports its own battery voltage, and we monitor to ensure that battery voltages are largely the same across 
devices throughout the experiment. 

B. Protocol 

1) Software: A NesC/TinyOS embedded program is written to operate the following protocol: 

Frequency Hopping: From the 902-928 MHz band, 14 center frequencies are chosen. The nodes are programmed such that 
each node hops across all the 14 frequencies in each cycle. The time duration between frequency hops is three seconds. 

Synchronization: Synchronization is required so that frequency-hopping sensors are all transmitting and receiving on the 
same frequency at the same time. One of the frequencies in the frequency band is considered as a synchronization frequency 
and is repeated three times each cycle so that neighbors can synchronize with each other more quickly. The dwell time on 
each frequency includes a period in which all sensors transmit a packet and receive packets from other sensors, and a period for 
switching frequency. 

Pairwise measurements: Each node measures path loss on all links with all other nodes at each frequency. A TDMA-based 
MAC scheme is used in which each node broadcasts its pairwise measurements during an assigned slot, to avoid interference. The 
data sent by a node in its packet transmission includes the RSS values recorded during the previous period, a unique sequence 
number, its transmit power, and its battery voltage. 

2) Receiver Base: The receiver base is a mica2 node connected to a laptop, loaded with a special receiver program which 
synchronizes to the frequency hopping schedule of the nodes and communicates all the received packets serially to the laptop for 
storage and later analysis. 

IV. Experiment 

This section describes the use of the NCMS described in Section |III] to measure a network deployed in an ensemble of 15 
different environments. These measurements will enable the statistical analysis and model development in subsequent sections. 

A. Motivation 

Ideally, statistical characterization of the radio channel for multi-hop networks would proceed as follows; deploy K networks, 
each with A^ nodes positioned with the identical geometry in the same type of environment, but each network in a different place. 
For example, one might deploy the NCMS in a grid, in K different office buildings. 

In reality, its not economical to carry out the measurement campaign in K different office buildings, mainly because it is 
difficult to obtain access to carry out measurement in many different office areas, and it is difficult to position sensors in exactly 
the same geometry without moving obstructions to make space for each node. If the environment must be altered to measure it, 
we might as well randomly alter the entire environment. 



In fact, in this campaign, we occupy a single environment and randomly vary the object locations in that environment. We 
start by deploying nodes in an empty classroom in the Merrill Engineering Building at the University of Utah. A 4x4 square 
grid of mica2 nodes is set up with 4 ft (1.22 m) separation between neighboring sensors. Within this deployment area, different 
arrangements of obstructions are randomly generated. 

1) Random Environment Generation: For reasons of portability, the obstructions used in this campaign were cardboard boxes 
of size 61 cm x 41 cm x 61 cm (24 in x 20 in x 25 in). In order to make the boxes significant RF scatterers, we wrap the 
cardboard boxes with aluminium foil. Foil-wrapped cardboard boxes represent metal obstacles which might be present in office 
environments. 

We generate (in Matlab) random positions for 10 boxes to be placed in the area of the deployment. The Matlab script is written 
to ensure that boxes do not lay on top of any of the 16 sensors (which are placed on the floor). Beyond that restriction, the 
rectangular boxes may be positioned anywhere in the environment and may be positioned with either N-S or E-W orientation i.e., 
with their longer sides parallel or perpendicular to the X-axis as shown in Fig. [3] 

B. Experiment Procedure 

After random placement of the 10 obstructions, the campaign proceeds by powering on the 16 nodes and receiving and recording 
the measured path loss data in a file on a laptop. Each node runs the algorithm described in Section |III] After 10 minutes of run- 
time, the nodes are turned off. The process continues with the next measured network by randomly changing the obstruction 
locations and repeating the experiment. Fifteen network realizations are measured in this manner. 

V. Statistical Analysis 

This section presents the statistical analysis of the data collected by campaign described in Section |IV] We first estimate the 
path loss model parameters of© and (|2|i. Next, we analyze the shadowing loss correlations which exist on different pairs of links. 

A. Analysis of Received Power 

We denote the number of the deployment experiment as to S {!,..., Af }, where M is the number of deployments (here, 
M = 15). We denote the set of frequencies measured as ^. The received signal power between node i and node j for experiment 
771 at center frequency / G S^ is denoted i^ " (/) and can be written using Q and Q as 

P^f (/) = Pt, - no - lOn, log Ij - X[f - Y}f if) , (4) 

where Y^"^' (f) is the non-shadow fading and X>"^' is the shadow fading on link {i,j) during experiment m. Shadow fading is 
considered to be constant across the frequency band, as discussed in Section |T] We denote the frequency average received power 

p{m) ^ J_ V^ P^"^h f\ 



From (|4ll, we can write P ^ '^ as 



Pij) = Pr, - n„ - lOn, log ^ - 47) -^Y. Y'f (/) ■ (5) 



In other words, (|5]l can be written as. 



P^ = Pt, - Ho - lOn, log |j - X^ - Y^ , (6) 



where X; " = m X) f es ^ij (/)■ Because Fj " is an average of measurements at many different frequencies, we argue that it 
may be well-represented as Gaussian (in dB), regardless of the underlying frequency-selective fading mechanism {e.g., Rayleigh 
or Rician). Since X^^' is also log-normal 1201 . we expect the sum Z^^' to also be Gaussian (in dB). 

A linear regression of the frequency averaged received signal powers {pI™^ }i,j versus known distances {dij}ij is used to 
estimate the constants [Pt — Ho) and rip ^ for each experiment m. In our experiments, we have used Aq = Im. Since all nodes 
are set to the same transmit power and have approximately equal battery voltages, and since we estimate {Pt — XIq) in addition 
to Up, we are not required to know the exact transmit power Pt at the current battery voltage of the nodes in the network during 

experiment m. The linear regression also determines the variance of Z^^. 



B. Analysis of Link Correlations 

In this subsection, we describe the computation of the correlation in fading between pairs of Hnks. This requires computing 
correlation in the sample values of Z^"^ for different pairs of links {i,j) as described in Section [V-AI 

1) Similar Geometry Links: We use the term "link geometry" to describe for two links, link a and link b, the relative coordinates 
of the end points of the two links. In a grid network, there can be many pairs of links with the same link geometry (within a 
rotation). As one example, the link pair of link a and link b shown in Fig.|2] is repeated 16 times in the network as shown in Fig.|4l 

Let L denote the number of times a particular link geometry is repeated in the network. We denote the p link pair as the two 
links (ip, jp) and (fcp, Ip), where p <E {1, . . . , L}. Then 2'^ and ^j.™; , where m G {!,..., M}, represent the total fading on 
the p repeated link pair for experiment number m. Then vectors Zi™ and Z^™ are defined as 

We then define vectors Z^ and Z;, as 

Z, = [ZW^...,Z(^^)T , Z, = [zW^...,ZrY. (8) 

Vectors Z^ and Zt are both LAIxl sized vectors. Together they contain all measured total fading values for pairs of links which 
share a particular link geometry. The correlation coefficient of total fading on link a and link b, pZa,Zb > can be computed by taking 
vectors Z^ and Z(, as sample values of total fading for link a and link b respectively. 

We have computed the correlation coefficient for total fading on link a and link b for a variety of link geometries. TableUshows 
the results for various link pair geometries. We also run a hypothesis test to determine if the measured correlation is statistically 
significant. This test compares hypotheses. 

Ho : Za and Zb have p = 0, 
Hi : Za and Zf, have p ^ 0- 

We report P [measuring p\Ho] using the method described in f2ll pp. 427-431], in Table|T] The proposed correlated link shadow- 
ing model and the Gudmundson model, also mentioned in TablelH will be discussed in SectionlVTl 

C. Discussion 

The results show that, for many link pair geometries, it is extremely unlikely that the fading losses measured on the pair of links 
are independent. For 15 of the 28 studied link geometries, there is statistically significant non-zero correlation. Those 15 links are 
consistently those geometries in which the two links are proximate, i.e., their lines from transmitter to receiver partially overlap, 
or nearly overlap. The likelihood that the measured correlation coefficient was measured by chance in the case when p = is 
extremely small, i.e., less than 0.5%, for 11 of the 15 link geometries which showed correlation. 

Also note that the correlation coefficients are relatively large in magnitude. The highest p is 0.33, six link geometries have 
p > 0.20, and eleven link geometries have p > 0.10. Fading loss on one link is obviously not purely determined by the losses 
experienced on its geographically proximate links; however, the correlation coefficient indicates that knowing the losses on the 
proximate links can give quite a bit of information about the loss on that one link. 

VI. Joint Path Loss Model 

In this section we present a model to describe the experimentally observed characteristic of correlated hnk shadowing. We start 
with the assumption that shadowing loss experienced on the links in a network is a result of an underlying spatial loss field p(x), 
such that shadowing on a link is increased when its path crosses areas of high loss p(x). We show how this assumption results in 
agreement with existing path loss models when considering a single link. We then show how it accurately represents correlated 
shadowing losses when jointly considering links in a multi-hop network. 

A. Shadow Fading Model 

In particular, we assume that the underlying spatial loss field p(x) is an isotropic wide-sense stationary Gaussian random field 
with zero mean and exponentially-decaying spatial correlation. The covariance between p at positions xi and X2 as 

£;[p(xi)p(x2)] =i?p(xi,X2) = i?p(||x2 -xill) = ^cxp r^^^^M") . (9) 

where ||x2 — xi|| is the Euclidian distance between xi andx2, S is a space constant and ax is the standard deviation of the shadow 
fading. The contour plot of a realization of such a random process is shown in Figure|5] 



Many mathematically valid spatial covariance functions are possible ll22ll . We justify the use the covariance function in (|9]l 
because of its basis in a Poisson spatial random process. Poisson processes are commonly used for modeling the distribution of 
randomly arranged points in space, and we suppose that attenuating obstructions might be modeled in such a fashion as well. 
Without detailing a specific model for the spatial extent or value of attenuation of each obstruction, we note that many Poisson 
processes (or derivatives of Poisson processes) have covariance functions with an exponential decay as a function of distance, as 
®. 

We propose to model the shadowing on all links as functions of the spatial loss field. We propose to model the shadowing on 
link (m, n), X„i^n, as 

Xm,n - T ^^77 / " P(x)dx. (10) 



X^i X^i 



1/2 



a) Single-Link Properties: This model agrees with two important empirically-observed link shadowing properties: 
Prop-I The variance of dB shadowing on a link is approximately constant with the path length 1201 . 131. 141. 
Prop-II Shadow fading losses are Gaussian. 
The model in (fTOl i can be seen to have Prop-II, since Xm,n is a scaled integral of a Gaussian random process. 
The proposed model has Prop-I when ||xj — Xi|| >> 5. We show this by considering the variance of Xa, 

Var [Xa] = E[Xl] 

1 n r' ^ n,^ n^ , T,. (H) 



/ / Rp{\\(3-a\\)da'^d(3. 



(12) 



Using (|9]l as the model for spatial covariance, (fTTT i is given by 

Var [Xa] =a\\l+ „ ^_^ e-"-^— 'H/^ 

When ||xj - Xj|| >> 5, 

Var [Xa] « a\. (13) 

b) Joint Link Properties: Next, consider two links a = (i, j) and h = (fc, /), as shown in Fig.[5]with shadowing Xa and X^, 
respectively. Consider the covariance of Xa and X^, 

Co^{Xa,X,) = - -^ -— j j e-'^^^doFdl3. (14) 

where C„i,„ is the line between points x„i and x„. Since £'[X(j] = -B[^&] = 0, the correlation coefficient between Xa and X^, 
Px^.x.As 

Cov{Xa,Xb) 

"''-''' v/Var[X„]Var[X,] ^^^^ 

^^-^^ "^ Allx- X-||l/2||xz X,l|l/2 / / ^ ' dcx d(3. 

0||Xj-x_,|| / ||XA;-Xi|| / Jc,jJCk,i 

The solution to ( fTSl ) is tedious to analytically derive. We use numerical integration to compute the value of px^,Xt^ and Matlab 
calculation code is available on the authors' web site Ii23il . 

B. Total Fading Model 

Since shadowing loss Xij is only one part of the total fading loss Zij ~ Xij + Yij, we must also consider the model for 
non-shadowing losses Yij . It is worthwhile to note that shadow fading and non-shadow fading are caused by different physical 
phenomenon, and thus Xij and Yij can be considered as independent. The variance Var [Zij] is thus 

a^B - Var [Z,j] = Var [X,,, + Y,J - Var [X,,,] + Var [Y,,,] . (16) 

Non-shadow fading is predominantly composed of frequency-selective or small-scale fading, which can be well-approximated to 
have zero correlation over distances greater than a few wavelengths. Since multi-hop networks typically have sensors spaced more 
than a few wavelengths apart, {i^^.j} are considered independent in this paper 

Thus the correlation coefficient between the total fading on links a and b, Za and Z5, is 

Coy {Za,Zb) 

PZa^Zi 



v/Var [Za] Var [Z^] 

r 1, ifa = 6 (17) 



^ V -^Px^,x, ~ -^Px^^x,, if a^b 



Equation ( fTTI ) indicates a linear relationship between the correlation coefficient of total fading and correlation coefficient of 
shadow fading. The correlation coefficient, pz^.z^,, is the measured p computed in Table J] The total fading variance a\Q was 
determined by the regression analysis in Section [V-A| 



C Estimation of model parameters from measurements 

Both the space constant 5 and the variance of shadowing a'jr must be determined experimentally from the data set. Specifically, 
we find the {S, a\) pair which best explains the correlations which exist in the link measurements. In other words, the goal is to 
find the value of [5, (j\) that results in highest agreement between measured and model-based correlation values. 

To accomplish this model fitting, we compute the model correlation, px^.x^, for a range of (5 e [0.1, 0.4] using ( fTsT i, for each 
of the 28 link geometries considered in Table I] At a particular value of 5, we compare the model correlation px^.x^ with the 
measured correlation pZa.z^ using linear regression. This linear regression returns a correlation coefficient, pc, which quantifies 
how well the model (using 6) agrees with the measurements. The highest value of pc uses 6*, the optimum 6 which matches the 
model to the measurements. Fig.|6]plots the correlation pc for 6 E [0.1 0.4]. We can observe that the curve attains the maximum 
Pc at S* = 0.21. The value of a^ is then determined from 6* using ([TtI i. and we see that crx/<^dB ^ 0.29. 

In summary, we have determined the two parameters of the correlation model, {S, ax) using the measurement data set. 

D. Comparison with Gudmundson Model 

In this section, we compare the proposed model of shadow fading correlation with an application of an existing model |fT6ll . 
Gudmundson's model addresses cellular radio networks where a mobile receiver (low antenna) communicates with a base station 
(high antenna). As the mobile receiver changes position with respect to the base station as shown in Fig.|7] there can be significant 
correlation in shadowing on the links with the base station. For a mobile receiver moving with a velocity v, and sampling signals 
at every T seconds, the correlation in shadowing Rx{k) is given as 

i?x(fc) = cri-a''"' where, a = ef^''. (18) 

where D is the reference distance, eu is the correlation in shadowing on Unks when the mobile receiver moves a distance D, and 
tj^ is the variance of the shadowing on a link. 

1) Application to Multi-hop Networks: Because the model of (fTsT l was not designed for ad hoc networks, it can only be applied 
to pairs of links which share a common node. This is a major limitation of the Gudmundson model which requires development 
of a new shadowing correlation model for multi-hop networks. Regardless, we consider here the application of (fTSI l to pairs of 
links which share a common node. The shadowing correlation between the two links from a common node to two nodes at x^ and 
Xj can be written as 

n r \ 2 l|xi— x,||/d /-ins 

i?x(x,,Xj) = axe'o ■ (19) 

Taking the logarithm of (fT9l l. we get a linear equation in (|x, — Xj ||, 

log Rx (x, , X, ) = log aj, + M^M log e„ . (20) 

The constants aj^ and e,, can be determined by running a linear regression between log Rx (x^ , Xj ) and measured correlation 
values (in Table|I|. 

Another limitation of applying the Gudmunson model to multi-hop networks is that it ignores the location of the common node. 
For example, in the two examples in Fig. [8] the correlation predicted by Gudmundson's model would be identical for both (a) 
and (b). Experimentally, the correlation varies significantly, from 0.21 in (a) to 0.05 in (b). Gudmundson's model is based on the 
assumption that the distance between the base station and mobile station is large compared to distance moved by the mobile. This 
assumption is not generally applicable to multi-hop networks. 

Table HH compares the ability of the proposed and Gudmundson's model to predict the measured correlation values. For the 
proposed model, we compare the p value from the 'Prop. Model' column of Table U with the 'Measured' column, for all 28 link 
geometries tested. For Gudmunson's model, we compare the 'Gud. model' with the 'Measured' p for the 21 link geometries to 
which the model can be applied. We observe that the measurements have 80.4% agreement with the proposed model, compared 
to 64.4% with Gudmunson's model. Note that while both models are 'fit' to the data, the comparison is valid since both models 
require fitting of two parameters {ax and S in proposed model and ax and e^ in the Gudmundson model) to the data. 

VII. Application of Joint Model 

In this section, we study the effects of shadow fading correlation in two fundamental multi-hop network examples, paths in three 
and four node ad-hoc networks. We show by analysis and simulation that the probability of a path failure can be significantly 
higher when links have correlated, as opposed to independent, link shadowing. 

To simplify the analysis we assume 



1) Packets are received if and only if the received power is greater than a threshold 7, and 

2) No packets are lost due to interference. 

These assumptions do not limit the results in this section. In fact, performance in interference is also affected by joint path losses, 
and is further impacted by correlated shadowing. 

We denote the normalized received power above the threshold for a link (r/i, n), as (3m,n, 

/5„.„ = :^ii^i^ (21) 

where 7 is the threshold received power and Pm.n is received power given in (|2]i. Link (m, n), by assumption, is connected if and 
only if P,n,n > 0. An important system parameter is the expected value of /3m, n^ 

Pm,n = E [Pm,n] ~ ' (22) 

(^dB 

where P{dm.n) is given in ([T]i- Intuitively, ^m,n is the number of standard deviations of link margin we have in link [m, n). If 
we design the multi-hop network with higher Pm,n, we will have a higher robustness to the actual fading in the environment of 
deployment. For example, one could set the inter-node distance to ensure that /3m „ = 2, and then link [m, n) would only be 
disconnected if total fading loss was two standard deviations more than its mean. 
We define two events relating to the connectedness of links, 

A = {Link (i, k) is connected} ~ {Pi^k > 0} 

B = {Link {i,j) and Unk {j, k) is connected} = {/?ij > 0} n {fij^k > 0}. 

Then AU B is the event that two nodes i and k can communicate, either directly or through an intermediate node j. We call the 
probability that node i and k cannot not communicate as the probability of path failure, 

1 - P[AU B] = I - [P[A] + P[B] - P[Ar\ B]]. (24) 

A. A Three Node Multi-Hop Path 

Consider the simple multi-hop path shown in Fig.|9ja), which represents a part of a typical multi-hop network. In this example, 
||xi — Xjll = ||xj — Xfc||. For node i to transmit information to node k, the message packet can take two routes. One is the 
direct link {i, k) and the other is a two hop path through a relay node j i.e., through link {i,j) and through link {j, k). If for our 
particular deployment, the link (z, k) fails due to high shadowing, there is a chance that the message can still arrive via links (i, j) 
and (j, k). This section shows that this 'link diversity' method is not as robust as would be predicted assuming independent link 
shadowing. 

From (I22I 1 and ([Til, the relationship between (3ij, (ij^k and /3i.fc is 



(ii,j = /3j,fe; and (ii^k = Aj - «• (25) 



lOWp logio 2 



where n 

According to the definition ( 1231 ), the probability of event A is, 

P[A]=P m,k > 0}] = Q {-P,x) = 1 - Q [hj - ^) (26) 

where Q (■) is the complementary CDF of a standard Normal random variable. 

1) Case ofi.i.d. Shadowing: Under the assumption that the shadowing across links in a network is i.i.d., the probability of 
event B is 

P[B]^P [{A.J > 0} n {(i,.k > 0}] = (1 - Q {h,)f. (27) 

From dZTb and ( l24b . the probability of path failure is 

l^P[A\JB]=Q (A,, - «) Q (A,, ) [2 - Q (A J ) ] • (28) 

2) Case of Correlated Shadowing: From the correlation values reported in Table HI we know that links (i, j) and (j, k) of 
Fig- IS a) are nearly uncorrected. Thus, the probability for event B is approximately the same as in i.i.d. case. The probability of 
path failure in this correlated case is derived in the appendix to be 



i-P[Ar\B] = i 



Q 



-A^i 



'l-Z'^.^.x., 



dA,fc, (29) 



where, 

Ml = A J + (A,fc - A J + i^)pxi,j,Xi, 



B. A Four Node Multi-Hop Network 

Next, consider the four node link shown in Fig. |9jb). For this linear deployment we assume ||xi—Xj|| = ||x_,— Xfc|| = ||xfc— Xi||. 
For node i to communicate with node I, the message packet can be routed in four ways as shown in Fig. IHb). An analytical 
expression for the probability of path failure is tedious, so instead we simulate the network shown in Fig. |9jb) in both the 
correlated and i.i.d. link shadowing. We take 10'' samples of the normalized received powers under both correlated shadowing 
and i.i.d shadowing models. We then determine from the result the probability that there is no path from node i to node I. 

C. Discussion 

We compare the probability of path failure between node i and node k for both the cases of i.i.d. and correlated link shadowing 
in Fig. [To] The analysis shows that when a multi-hop network is designed for Pij — 2, then the probability of path failure is 
120% greater in correlated shadowing as compared to i.i.d. shadowing. Increasing the reliability of the network by designing it for 
higher Pij only increases the disconnect between the two models. It is only when we design the network for very unreliable links 
{e.g., Pij = 0, for which link {i,j) is connected 1/2 the time) that the models have a similar result. Clearly, path connectivity is 
much more likely under the i.i.d. model than under the realistic correlated link shadowing model. 

The four-node example shows that as paths become longer, it becomes increasingly important to consider correlated link 
shadowing. While the 3-node network had a 120% increase in probability of path failure, the 4-node network showed a 200% 
increase in the same probability. While Figure [TO] show the results up to Pij — 2.5, higher values correspond to higher reUabiUty 
links, and reliable networks will be designed with even higher link margins. When networks are designed for high reliability, the 
effects of ignoring hnk correlations are dramatic. 

VIII. Conclusion 

A statistical joint path loss model for multi-hop (sensor, ad hoc, and mesh) networks is presented that relates the shadow 
fading on different links in a multi-hop network to the underlying shadowing field caused by an environment of deployment. A 
network channel measurement system is used to measure a multi-hop network deployed in an ensemble of environments. The 
data set is used to demonstrate and quantify statistically significant shadowing correlations among different geometries of links. 
The measured correlations agree with the proposed model, and can be applied to a greater variety of links than possible using 
an existing correlated shadowing model. Finally, this paper analyzes path connectivity in simple multi-hop networks to show the 
importance of the consideration of shadowing correlation when designing reliable networks. The probability of path failure is 
underestimated by a factor of two or higher by the current i.i.d. shadowing model. 

Future work will test other ensembles of deployments, both indoors and outdoors. The effects of correlated shadowing will have 
impact on higher layer networking protocols and algorithms, and in interference and multiple-access control, and future work will 
quantify this intuition. 

APPENDIX 
Here we present the derivation of the probability P[AC)B] in ( |29] l. From (ISTT i. we can note that Pij, jSj^k and l3i,k are 
joint Gaussian random variables. Thus the conditional distributions, f{Pi,j\Pi,k ~ b) and f{Pj^k\Pi,k = b), are Gaussian. The 
links {i,j) and (j, k) in Fig. |9]are observed to have very small or no correlation between them. Thus the joint distribution, 

fiPij, Pj.k\0i,k = b), can be approximated as: 

/(A-,,/?,-fe|A,fe = 6) « /(A,,|A,fc = &)/(/?,^fe|A,fc = b). (30) 

The joint distribution of /3ij, (3j^k and j3i^k is: 

f{|3^,J,pJ.k, (i^,k) = f{P^,0^,k = b) f {P,.k\[3^^k = &)/(ft,fc). (31) 

The probability P [^ n S] can be written in terms of joint distributions as: 

p[Ac^B] = P [{A,, > 0} n {/3,, ,. > 0} n {A,,. > o}] 

.f{P^M.k = b)J{Pj.k\l\k = b)f{p^,k)dP^,jdp,,kdp,^k 

-db. (32) 



{ft,j>0} J{l3j,k>0} J{l3i,k>0} 

2 (b-fii,k)^ 



lb>0 

where. 



/ [Q(-m/7l-Px. 
Jb>o ■- ^ ^ 



J -^i.k 



fii 4 E[{A,,|A,fc = b}] = A,, + ib- h, + ^)px.,,,x„,. 
The square in the RHS of (|32] | comes from the fact that for the link geometry considered, px k-Xi k ~ PXi -.Xt k- 
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Fig. 1 

Graphical depiction of (a) circular coverage model, and (c) coverage in the i.i.d. log-normal shadowing model, compared to the 
common depiction of (b) in which coverage area is a random shape. in (a) and (b), nodes are connected if and only if they are 
within the gray area, while in (c), nodes are connected with probability proportional to the shade (darker is more probable). 
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Fig. 2 

Example of factor in shadowing loss correlation. Because link a and link h cross the same environment, their shadowing losses 

TEND TO BE CORRELATED. 
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Fig. 3 

One realization of the random locations of the boxes among the 16 node locations (x). Each box (grey rectangles) occupies two 
pixels of this graph and can be placed either parallel or perpendicular to x-axis. 
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TABLE I 

Link Geometry and Correlation Coefficients (Observed, Proposed Model, and Model of [Gudmundson 1991]) 




Fig. 4 

Link pairs with identical link geometry in a grid deployment. Each link pair is shown with one link as a dotted (- -) line and 

another link as a solid lines (-). All link pairs with identical link geometry are shown. 





CoiTelation with Measured Data 


Proposed Model 


0.804 


Gudmundson's Model 


0.644 



TABLE II 

Comparison between the proposed model and Gudmundson's model 
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Fig. 5 

A LINK PAIR IN AN UNDERLYING SPATIAL LOSS FIELD 
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Fig. 6 
Variation of pc with <5. 
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Fig. 7 
Example of the motion of mobile receiver and base station position in Gudmundson's model 
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(b).' -^ -^ • 

Fig. 8 

A CASE OF TWO DIFFERENT TYPES OF LINKS, SHOWN BY (A) AND (B). THE GUDMUNDSON'S MODEL PREDICTS IDENTICAL CORRELATION FOR THE TWO 
CASES WHILE THE PROPOSED MODEL DOES NOT. EXPERIMENTALLY, THE CORRELATIONS VARY SIGNIFICANTLY FROM (A) 0.21 TO (B) 0.05 





(b) 
Fig. 9 

Example multi-hop networks of (a) three nodes and (b) four nodes. 
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Fig. 10 
Plot showing the variation of percentage increment in the P [link failure] for a 3 node and 4 node multi-hop network with 

normalized RECEIVED POWER OF THE SHORTEST LINK (i, j), ftj. 



